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Abstract 

Using a newly suggested algorithm of Gozzi, Renter, and Thacker 
for calculating the excited states of one dimensional systems, we de- 
termine approximately the eigenvalues and eigenfunctions of the an- 
harmonic oscillator, described by the Hamiltonian H = + gx'^. 
We use ground state post-gaussian trial wave functions of the form 
^(x) = A^exp[— blxp""], where n and b are continuous variational pa- 
rameters. This algorithm is based on the hierarchy of Hamiltonians 
related by supersymmetry (SUSY) and the factorization method. We 
find that our two parameter family of trial wave functions yields ex- 
cellent energy eigenvalues and wave functions for the first few levels 
of the anharmonic oscillator. 
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The central idea of the factorization method for ordinary differential equa- 
tions, introduced by Schrodinger [|l| and by Infeld and Hull to solve prob- 
lems in quantum mechanics, is the recognition that once the ground state 
energy and wave function of a one-dimensional potential problem are known, 
then the potential is determined, as well as the factors of the Hamiltonian. 
Recently, Gozzi, Renter, and Thacker proposed a simple strategy for solv- 
ing such one-dimensional potential problems approximately in a variational 
scheme, based on utilizing the hierarchy of Hamiltonians that are related by 
supersymmetry (SUSY) and factorization 0, H, |^. We apply their method 
here to the case of a one- dimensional anharmonic oscillator, using a simple 
variational wave function, and examine the results. 

First let us review briefly the formalism for describing the hierarchy of 
isospectral Hamiltonians. Consider the Hamiltonian Hi having k bound 
states whose discrete eigenvalues and eigenvectors are denoted by E^"* and 
'^^n' ■ Once we subtract the ground state energy, we can factorize the Hamil- 
tonian as follows: 



where the superpotential. 



Wi{x) = -j^[M'^\x)l (2) 

and 

V,^{l/2){W!-W[). (3) 
The second Hamiltonian H2 in the hierarchy, the SUSY partner of Hi, 

H, = ^AiA{ + E« (4) 

has the same spectrum as Hi, except that there is no state in H2 correspond- 
ing to the ground state of Hi. Introduce the notation that in E!;^\ n labels 
the energy level and (m) refers to the m-th Hamiltonian Hm- 
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One can show 0] [|| that the energy eigenvalues and eigenf unctions of 
the two Hamiltonians Hi and H2 are related by 

^Si = 2-^/^(4^1 - i?S^V^/^4*f (5) 

and in particular: 

+ H.^)vi>(2) = 2i/2(E{^) - (6) 

(JiJb 

so that the first excited state of the original Hamiltonian Hi can be ob- 
tained from the ground state of the supersymmetric partner potential. If 
we determine the ground state wave function \i/o and ground state energy 
{=E^i^) for H2, then we can repeat the above argument and determine a 
third Hamiltonian H^ as the SUSY partner of H2- Namely we can write 

H2 = ]^AiA\ + E'i'^ = ]^A\A2 + ^ (7) 

where 

A2 = ^ + 1^2(X). (8) 

Then 

H^ = ^A24 + Ei'^ ■ (9) 

We also notice that the Hamiltonian given by H2 = \AiAi has a ground 

state energy which is the energy difference e[^^ — Eq^\ Furthermore we have 
the relationships: 

Ei% = 41 = El^\ (10) 
^Si = 2-^2(4^^ - i5;«)-V24^(2) 

= (l/2)(4,i-i?^VV2(4l)^_^(l))-l/244^(3) (11) 

In particular, for = 1 we determine the second excited state of Hi from 
the ground state of H^. Clearly this method can be generalized to obtain k 
Hamiltonians of the hierarchy, where k is the total number of bound states 

of iJi. 
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The scheme of Gozzi et. al. [Q] is to determine the ground state wave func- 
tion \l/o^'' and energy e'q'^ of the original Hamiltonian Hi using a variational 
method, then to approximately find the superpotential 

W^i.(^) = -^{ln<^(x)}. (12) 

This allows one to approximately determine the first SUSY partner potential 
corresponding to H2 

V2.{x) = \[wl{x) + W[,{x)\ (13) 
obtained from the approximate partner Hamiltonian 

H2v = ^Ai,Al„ Ai, = ^ + Wi,. (14) 

whose ground state energy determines approximately the energy difference 
E^i^ — Eq^^ as discuused earlier. One next finds the variational ground state 
wave function and energy of the approximate SUSY partner Hamiltonian H2 v 
from which one can also calculate the approximate first excited state of the 
original Hamiltonian, and we can proceed as described above, using a new 
variational approximation for each ground state, until all the wave functions 
and energies of this approximate hierarchy of Hamiltonians are known. 

The interesting question is how accurate is such a scheme? All the Hamil- 
tonians of the variational hierarchy are not the real Hamiltonians of the hi- 
erarchy but approximate ones related to the variational approximation of 
the ground state wave function of the original problem. Thus we do not 
a priori know that we are calculating accurately the energy eigenvalues of 
the original Hamiltonian. We only know that we are accurately determining 
the energies of the approximate Hamiltonians of the hierarchy. Thus, for 
this scheme to work, one needs an extremely accurate method of finding the 
ground state wave function and energy of Hi. Also in this scheme the higher 
wave functions of Hi are obtained by repeated derivatives of the ground state 
wave functions of the approximate hamiltonians. Thus we expect a loss in 
accuracy as we go to the higher wave functions. 

In some recent work on soliton dynamics in various nonlinear systems 



including the nonlinear Schrodinger equation [10|, we have demonstrated the 
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efficacy of post-gaussian time dependent variational wave functions of the 
form 

u{x,t) = N{t) exp[-b{t)\x - . (15) 

These wave functions, when used to minimize the action, gave excellent global 
fits to various solitons as well as very accurate determinations of the soliton 
energies. We found that these wave functions arc globally accurate to a few 
percent or so, and give energy eigenvalues accurate to about 0.1%. 

Thus, we propose to use as our trial wave functions for the ground state 
wave functions of the hierarchy of Hamiltonians, 



Nk exp 



r 1 






-^2 


Pk 





1 



2nk' 



-1/2 



(16) 



We first scale the Hamiltonian for the anharmonic oscillator. 



(17) 



by letting x — > xjg^l^ and H — > g^/^H. Then we find the ground state 
energy of the anharmonic oscillator and the variational parameters pi and 
ni by forming the functional 



1 (f 

£;o(pi,ni) = (*o|-2^ + ^1*o: 



(18) 



We determine pi and ni by requiring 

0, 



0. 



dpi drii 
The equation for the energy functional for the anharmonic oscillator is: 

W (r. ^\ ^? ^(^~ 2^) , 2^(2^) 



(19) 



2pi 



V2ni- 



V 2ni - 



(20) 



Minimizing this expression, we obtain the variational result 

Eq = 0.66933, rii = 1.18346, pi = 0.666721 . 



(21) 
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This ground state energy is to be compared with a numerical evaluation 
which yields 0.667986, which is of the accuracy obtained in our soliton cal- 
culations. Since the trial wave function for all ground states is given by (]TE|), 
the variational superpotential for all k is 



W, 



k V 



2nfc 



(22) 



Since we are interested in the energy differences Ek — Ek~i of the anharmonic 
oscillator we will concentrate on the variational Hamiltonian 



Hk+lv — T^-^kv-^kv 



(23) 



which approximately determines these energy differences as discussed above. 

First, let us compare the exact potential — Eq with the variational one 
obtained using the ground state variational wave function: 



VUx) = -WUx)-W[,{x) 



(24) 



This is shown in the lower two curves of Fig. 1 for the most relevant range 
< a; < 1.5 for determining the ground state wave function. We see that 
the variational potential is quite different from the exact one, although it has 
almost the same ground state energy. It is also interesting to compare our 
variational "partner potential," 



(25) 



with the exact one, obtained by first solving the Riccati equation for Wi 

1 



Wl{x) - W[{x) 



X 



E, 



, 



or equivalently solving the Schrodinger equation for 

1 rf2 



2dx^ 



+ x^ 



^0 = ^0^0 



(26) 



(27) 



for ^^0; and determing Wi from the logarithmic derivative of ^E'o- The result 
of this calculation of the exact partner potential and its variational approxi- 
mation are shown in the upper two curves of Fig. 1. 
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We obtain the approximate energy splittings by minimizing the energy 
functional 



AE,{p„ n,) = 1( + Wl + WLl^l^'^ ) = - • (28) 

These integrals are trivial to carry out, and one obtains the simple recursion 
relation 



AEk{pk,nk) 



r (2 - nl^^ ( pk y^-' r 



r(^) 2p, \p,_,J r(^) 



2nfc 



+ ^(2n,_i-l)(^) ^ ;\\ ^ . (29) 

2p. \p,-.l r(4) 

One can perform the minimization in p analytically leaving one minimization 
to perform numerically. 

The results for the variational parameters and for the energy differences 
are presented in Table 1 for the first five energy eigenvalues and compared 
with our numerical calculation, based on a shooting method, which essentially 



agrees with the results of Hioe and Montroll|]TT|, who made an exhaustive 
study of the anharmonic oscillator. For the ground state and first two excited 
states, we get excellent agreement with the numerical results. However, for 
the third excited state and higher, the variational approach relaxes to a 
Gaussian which leads to equal spacing with a AE of around 2.3. Thus 
for the higher states our method fails to give the correct large n behavior, 
which from WKB results, is known to be 1.376(?t, + \Y^^- So this method 
works where WKB fails and vice-versa. Of course, our results for the higher 
energy levels would be improved by using more variational parameters, which 
however would spoil the simplicity of the calculation presented here. 

Finally, we have calculated the approximate excited states by repeated use 
of the approximate node insertion operators. The results for the first three 
variational wave functions are compared to the exact numerical results in Fig 
2. We see that we are doing an excellent job in the regime < x < 1.5, which 
is the most important region for determining the energy eigenvalues. The 
higher eigenvalues require more and more derivatives, and thus the accuracy 
of the results diminish for large x where the variational potential is not so 
accurate. 
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Table 1: Variational Parameters and Energies 



level 


n 


P 


En - 


En-l 


variation 


exact 





1.183458 


0.666721 


0.669330 


0.667986 


1 


0.995834 


0.429829 


1.727582 


1.725658 


2 


1.000596 


0.435604 


2.316410 


2.303151 


3 


0.999917 


0.434779 


2.297082 


2.638935 


4 


1.000012 


0.434894 


2.299820 


2.908578 
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FIGURE CAPTIONS 



Fig. 1: Exact (solid line) and variational (dashed line) potentials for the 
anharmonic oscillator (lower curves) and its SUSY partner (upper curves). 

Fig. 2: Exact (solid line) and variational (dashed line) wave functions for 
the lowest three levels of the anharmonic oscillator. 
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